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\^ ' Abstract 



Numerical investigation of the static spherically symmetric vacuum so- 
lution of the Logunov equations confirms the analytical results and demon- 
^~>) ' strates a strong repulsion at sub-Planckian distance from the Schwarz- 

^■f^ ' schild-like singularity, which is stable against the small perturbations. 

. We suppose that the final stage of the stellar collapse in the relativistic 

theory of gravitation can entail the gravitational bursts with the gravi- 
tons overproduction due to the strong deceleration in the vicinity of the 
Schwarzschild sphere. The rigorous investigation of this topic requires a 
O ' quantum-field consideration. 

?^ ■ 1 Introduction 

, . The relativistic theory of gravitation (RTG) is built as a classical field theory 

rS I in the Minkowski spacetime ^ |2 El 0]. Since the gravitational field is de- 

■ scribed by the symmetric second-rank tensor, it manifests itself as an effective 

Riemannian spacetime. Thus, the RTG is the bi-metric theory, which needs 
the gauge-symmetry violation to be physically self-consistent. This requires a 
massive graviton. As a result, the field equations contain both metrics (i.e. 
both Riemannian and Minkowski) and gauge-fixing condition, which eliminates 
the undesirable spin-states of graviton. The cosmological consequences of the 
non-zero graviton mass impose a limitation on its value, which cannot exceed 
10~66 g j^, but virtually can be incomparably smaller than this value 

The non-zero mass of graviton causes the repulsion (antigravitation) at small 
distances and the deceleration of the cosmological expansion at its late stage. 
For the black-hole physics is most important the antigravitation prohibiting an 
unlimited collapse of the massive stellar objects. This topic had been consid- 
ered earlier in the pure analytical framework (see, for example, \7\). However, 
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the RTG field equations (Logunov field equations) don't allow an exact ana- 
lytical solution even for the static spherically symmetric field and need some 
approximations in order to be solved in a closed form. Moreover, the physical 
consistency of the obtained solutions needs their stability at least against the 
small perturbations. 

In this work we consider the numerical solutions of the Logunov field equa- 
tions in the vicinity of the Schwarzschild sphere and compare their with the 
analytical approximation. This comparison demonstrates a high accuracy of 
the latter. Also, we analyze the stability of the analytical metric in the vicinity 
of the Schwarzschild sphere. Although this metric provides a strong repulsion 
for the radially falling particle, it is stable against the small perturbations. 

2 Numerical approach to the Schwarzschild-hke 
singularity 

The Logunov equations for the massive gravitational field have the following 
form [HlllElEI: 

- ^ (<5- + - ^C.^'^Ska) = -^r-, (1) 

D^r" = 0. 

Here g^^ is the metric tensor of the effective Riemannian spacetime, g'^'^ = 
\/—gg'^'^] 7'^^ is the metric tensor of the background Minkowski spacetime, 
is the Einstein tensor, Tj^ is the matter energy-momentum tensor, is the 
covariant derivative in the Minkowski spacetime; w? = (mgc/tif' , nig is the 
graviton mass; k is the Newtonian gravitational constant. 

Spherically symmetric static interval of the effective Riemannian spacetime 
can be written in the following form: 

ds^ = U (r) dt^ - V (r) dr^ - W {rf {d0^ + sin^ 9 dip^) . (2) 

Henceforward it is convenient to use the dimensionlcss quantities and we 
shall normalize the geometric units to the Schwarzschild radius = 2>cM / c?: 
W — X Ts and r = y r^- Then the field equations written in the interval 1^)1 

allow the system of the master equations — \ ( ^^^e™" ) ■, C — ^—^-^) Mi'- 
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= Cpx^ (3) 

Here we supposed that Tq = p, T™ = —p/c^ {m — 1..3). 
Using the new functions we write 




Since e is an extremely small quantity (see below), we can suppose 

(5) 

Then in the vacuum from (|3I4I5|I we have 
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Several approximations based on the assumption x ~ xi (where xi > 1 is 
the Schwarzschild-like singularity point) allow the analytical expression for the 
interval in the vicinity of the Schwarzschild sphere |3 IH] : 
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f + sin^ e dip') 



(7) 



The obvious differences from the Schwarzschild metric result in the main phys- 
ical consequence: the radially falling particle is decelerated in the vicinity of Xi 
and cannot reach the singularity. This statement is in accord with the causality 
principle in the RTG, which requires xq > x\ [xq is the position of the coUapsar 
surface) . 

To avoid the approximations required for the analytical expression ((JJ, sys- 
tem © was considered numerically. Since the upper threshold for the graviton 
mass is 1.6 x 10~^^ g 0, e has an extremely small value for the stellar objects 
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Figure 1: Dependence of the dimensionless metric functions on the dimensionless 
distance from the Schwarzschild sphere. Sohd - numerical simulation, dotted - 
metric Q , dashed - Schwarzschild metric with singularity at , dashed-dotted 
- Schwarzschild metric with singularity at xi. Point in (a) corresponds to 
analytical metric 



(e.g. e ~ 10"'"' for M = lOM©). This allows the assumption ^ = 1 and Eqs. 
(gEll result m: 

Far off xi or for e = Eq. (|SJ leads to the Schwarzschild metric. The 
behavior of U and V in the vicinity of xi for nonzero e is shown in Fig. 1 
(dashed curves present the Schwarzschild solution). The point in Fig. 1, a and 
the dotted curve in Fig. 1, h correspond to the analytical interval {T)) but for 
the numerically obtained value of xi. 

Firstly, we have to recognize a good accuracy of the analytical estimation 
in the vicinity of xi. Note also, that the Schwarzschild solution for V (dash- 
dotted curve in Fig. 1, 6) is quite accurate if = xi. Secondly, the differences 
between Schwarzschild and Logunov metrics appear at the extremely small dis- 
tance from (sub-Planckian for the stellar objects). It is clear that the latter 
results from the small value of rUg and requires the quantum consideration of 
the Schwarzschild sphere vicinity. The dependance of xi on e is shown in Fig. 
2. One can see that the deviation of xi from is very small even for the 
extra-massive coUapsars. 

The approach to xi of the free-falling particle is described in the vicinity of 
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Figure 2: Distance of the singularity point from the Schwarzschild sphere de- 
pending on the dimensionless mass parameter e. 
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Figure 3: Velocity of the falling particle depending on the distance from the 
Schwarzschild sphere. Solid - numerical result from the Logunov equations, 
dashed - result from the Schwarzschild metric. 
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the Schwarzschild sphere by the approximate expression [7]: 



dW (s) 
ds 



= -2 




X — Xi 



(9) 



which demonstrates a strong deceleration due to antigravitation induced by the 
massive graviton. Fig. 3 shows the numerical result obtained from Eqs. (|4I8|I 
and expression 



One can see, that the velocity of the falling particle (for the remote observer) 
approaches to the velocity of light (c = 1 in our normalization). However, 
in immediate proximity to xi the particle is decelerated by the short-range 
antigravitation and its velocity decreases up to zero (i.e. the Schwarzschild 
sphere is physically unreachable). Surely, as this deceleration takes a place at 
the extremely small distance, we need a quantum consideration. Nevertheless, 
one can suppose that the final stage of collapse can entail the over-production 
of the gravitons (some kind of the gravitational burst). Note also, that in spite 
of xi — 1 ^ 1 the time of the light propagation from r = lOOrg to the point, 
where the deviation from the Schwarzschild metric begins, is ^ 0.02 s for M = 
IOMq. Thus, the processes induced by the massive graviton in the vicinity of 
the Schwarzschild sphere are not hidden from us by the cosmological-scale time 
and we face the challenge of the stability of the observable massive coUapsar-like 
objects [in] . 

3 Stability of the spherically symmetric metric 

Although the problem of the coUapsar stability can not be reduced to the sta- 
bility of the metric against the small perturbation, the latter is its necessary 
condition. We shall base our analysis on the classical method of Ref. ^J. Q 
is chosen as the perturbed metric. 

The general perturbed metric has a form: 




(10) 



for the modified interval 



ds^ = U {W) dt^ 



V* (W) dW^ - {de^ + sin^ 6 dip^) . 



(11) 
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U{r), 0, 0, 

0, -F(r), 0, 

0, 0, -W{r)\ 

S uj {t, r,9)x S q2 {t, r,6)x S qs {t, r, 6) x 

W{rfsin^e ' W {rf d ' W {rf sin^ 9 



g^u = (12) 
5 UJ {t, r, 0) X 



W {rY sin^ 0, 
6 q2 (t, r, 0) X 
W {rf sin^ 0, 
5 q3 (t, r, 0) X 
W {rf sin^ 0, 

-W{rfsm^0 



where 5 is the smaU expansion parameter; uj {t, r, 0), q2 (t, r, 0) and (73 (i, r, 0) are 
the perturbation functions. 

Taking into account only hnear on 5 terms, using dimensionless quantities 
and supposing that xi — and 



UJ {t, r,0)=uj {r, 0) exp {-iat) , (13) 

92 (i, r, 0) = q2 (r, 0) exp (-icrt) , 

93 {t, r, 0) = 93 {r, 0) exp (-icri) , 

we can obtain from {T)) and second field equation 



3cos0q3{y,0) dq2{y,0) , dq3{y,0) 
7^2 (?/, 6') + 8(72 (y, 6*) y + 2y (y - 1) 



smf 



a?/ 



96* 



(14) 



2ai 



y^Lu {y,0) = O. 



Eq. H14() requires w (y, ^?) = if (72 {y,9), 93 (2/:^) are real and a is arbitrary 
and real. The latter requirement corresponds to the stability of the perturbed 
metric (|7I12|I in the framework of the linear stability analysis. More general 
consideration will be given elsewhere. 

Now we consider first field Eq. ^ for the perturbed metric. The obtained 
expressions are too awkward to be written here and can be found in [Sj. The 
axial perturbations are defined by the perturbed (1, 3) and (2, 3) - components 
of the field equations. Introducing the new function 



dq3iy,0) dq2{y,0) 



dy 



80 



sin'' 6', 



(15) 



we can obtain after some manipulations (see [S]) the equations defining the 
perturbation functions: 
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(16) 



= 293(2/,^) 
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e sm 



where A (y) = y^ (l — - ) . Now it is possible to eliminate (72 and ^3 . This 



results in 



d_ 
dy 



'l^[A{y)Q{y,efy 
Q{y,0) y sinHB(y) 



d_ 

de 



sin^ e A (y) 



Q{y,o) 



sm 



(17) 



Here A(y) = -2 + ^ - B (y) = -^wlzl and we used the smallness of e in 

order to eliminate tlie angular-dependence from A and B. 

Eq. H17(l permits the separation of variables: Q{y,0) = E{y)'i'{9). Then 
for the angular part we have the Gegenbauer's equation: 



cP COS0 d 



^ [9) + 2n^' {9) = 0, 



(18) 



where n is some number. Eq. (|18l) can be solved through Legendre functions 
and then 



a + 2) (^-1) 

n — , 

2 

I is integer. 

Taking into account that e <C 1, the equation for the radial part is 



(19) 



2eyA (y) -^E (y) + ey (4y - 1) (y) + (20) 

2S (y) [a^yS _ - 2ey + 3e] = 0. 

In immediate proximity to the Schwarzschild sphere and outside the S (y)- 
extremum (which appears as a result of the angular number growth, see below 
Eq. 12Hl 'l. we can neglect the first term in Eq. H2()(l . Hence, the approximate 
solution of Eq. ((201 is: 



S (A) « C (1 + eXf (3 + 4,A)(--^+i«^("-i°))/32. exp L^!ii±f^H^±M 

(21) 

where C is the constant of integration and A > 0: y = 1 + eA. 

Eq. 1)21(1 suggests that the appropriate asymptotic for S is provided only 
by > 0. Thus, the perturbation increment a is real in the agreement with 
the stability requirement. The stability far from xi, where the metric coincides 
with the Schwarzschild one, is provided by the stability of the latter. 
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4 Conclusion 



The numerical simulations based on the Logunov field equations demonstrate a 
high accuracy of the approximate analytical solution for the static spherically 
symmetric metric in the RTG. It is shown, that there exists a strong repulsion 
induced by the massive graviton at the Planckian distance. As a result of 
this repulsion, the collapse of the stellar objects has to entail the burst of the 
gravitational radiation. However, the rigorous analysis of the collapse requires 
taking into account the quantum effects. Some case for the stability of the 
coUapsar is provided by the stability of the static spherically symmetric metric. 

As the outlooks for the future analysis we can point out the following topics: 
i) more general stability analysis including higher-order perturbations; ii) classi- 
cal consideration of the possible static configurations of the matter in the static 
spherically symmetric metric; iii) realistic collapse models taking into account 
the energy loss due to gravitational radiation in the vicinity of the Schwarzschild 
sphere. 
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